Free Magnetic Moments in Disordered Metals 



O 
O 



oo 



Stefan Kettemann 1,2 and Eduardo R. Mucciolo 3 

1 Institut fur Theoretische Physik, Universitdt Hamburg, Jungiusstrafie 9, 20355 Hamburg, Germany, 

2 Max-Planck Institute for Physics of Complex Systems, Nothnitzer Strafie 38, Dresden, Germany, and 

3 Department of Physics, University of Central Florida, Box 162385, Orlando, FL 32816-2385, USA, 

(Dated: November 23, 2005) 

The screening of magnetic moments in metals, the Kondo effect, is found to be quenched with a finite prob- 
ability in the presence of nonmagnetic disorder. Numerical results for a disordered electron system show that 
the distribution of Kondo temperatures deviates strongly from the result expected from random matrix theory. 
A pronounced second peak emerges for small Kondo temperatures, showing that the probability that magnetic 
moments remain unscreened at low temperatures increases with disorder. Analytical calculations, taking into 
account correlations between eigenfunction intensities yield a finite width for the distribution in the thermody- 
namic limit. Experimental consequences for disordered mesoscopic metals are discussed. 
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In a metal with antiferromagnetic exchange interaction be- 
tween a local magnetic moment and the conduction electrons, 
correlations cause a change in the Fermi liquid ground state. 
The screening of the localized spin by the formation of a 
Kondo singlet below the Kondo temperature, Tk, is affected 
by disorder in various ways. Fluctuations in the exchange 
coupling due to random positioning of magnetic moments re- 
sults in a dispersion of Tk fH- Since Tk is defined by an 
integral equation similar to the BCS equation for the critical 
temperature of a superconductor, one could expect, by anal- 
ogy, the Anderson theorem |21 to be valid. In that case, the 
leading correction to Tk would be of order T^/g, where 

T^ is the bulk Kondo temperature in the clean limit and 
g = EpT, the dimensionless conductance, is the ratio between 
the Fermi energy Ep and the elastic scattering rate 1/r. How- 
ever, the equation for Tk requires a sum over the local density 
of states (LDOS) rather than the global one (DOS). It is known 
that disorder induces correlations in energy in the LDOS jj]- 
Thus, impurities cause additional fluctuations on the LDOS 
and thereby on Tk- 

Let us begin by describing the relevant approaches to the 
different regimes in the problem. Consider a single magnetic 
impurity in a closed, phase-coherent metallic grain where the 
energy levels are discrete, a Kondo box |4, 5], with mean level 
spacing A = l/vL d and Thouless energy E c = D e /L 2 . 
Here, L is the system linear size, v the DOS, D e the diffu- 
sion constant, and d the space dimension. We note that Tk 
does not exceed 1/r in metals, ruling out perturbation theory 
in 1/r. However, the regime E c < Tk < 1/r applies to 
a wide range of metals and an expansion in terms of diffu- 
sion diagrams is permitted in that case |@]. For smaller grains, 
when A < Tjc < E c , random matrix theory (RMT) yields a 
distribution of Tk that scales with A alone. When the grain 
is so large that the localization length £ is smaller than L, E c 
and A are irrelevant and Tk is determined by the average en- 
ergy level spacing in the vicinity of the magnetic impurity, 
A e = \/v£, d |0. Thus, one expects RMT to be applicable 
on the scale A^ in this local random matrix theory (LRMT) 
regime. However, corrections to RMT and LRMT due to cor- 
relations of wave functions of order 1 jg are found to be es- 



sential for determining Tk- Finally, when Tk is smaller than 
A or Ag, the distribution of Tk is mainly determined by the 
coupling of the magnetic impurity to a two-level system. 

In the symmetric Anderson model of a magnetic impu- 
rity 1 8], the exchange interaction matrix elements are esti- 
mated as Jki — At* k tu/U, where U is the on-site Coulomb 
energy and the hopping matrix elements connecting the 
localized state (f>i(r) with a delocalized state ipk( r )'- Uk = 
/ d d r 4>*(r) T^fc(r), where T is the kinetic energy opera- 
tor. For states localized at r = with a radius clq, Jki ~ 
J tpl(0) ipi(0) , where J = 4(l/m*a§) 2 /C/ is the exchange 
coupling, to* the band mass. In the antiferromagnetic case 
(J < 0), J is strongly renormalized. The temperature at 
which the second-order correction is equal to the bare value 
defines T K |9], 
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where N is the number of states in the band of width D. Here, 
we have used the rescalings x — D/J and n = 7k/ A, with 
x n = L d \ip n (r)\ 2 equals the probability density of the n-th 
eigenstate and s n = (E n — Ep )/ A is the eigenenergy relative 
to the Fermi energy in units of A. The number of electrons in 
the Kondo box is taken to be even. Equation Q is valid in the 
weak coupling regime, when Tk <C D. It coincides with the 
self-consistent solution of the one-loop renormalization group 
(RG) equation |8, 10] up to the tanh factor, which accounts 
for the finite-temperature occupation. While the approxima- 
tions involved in deriving Eq. Q fail to describe the Kondo 
system below Tk, Eq. Q does yield the correct low-energy 
scale in the effective Fermi liquid theory, where Tk deter- 
mines the Landau parameters [11]. Nevertheless, for finite 
systems, there may be deviations from one-parameter scaling 
due to sample-to-sample mesoscopic fluctuations 1 12]. 

In order to gain insight on the statistical properties of Tk 
over a wide range of values of J, we solved Eq. Q numer- 
ically and compared the resulting distribution and moments 
of Tk with approximate analytical expressions appropriate to 
the different regimes already mentioned. 
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FIG. 1 : (Color on-line) The distribution of Tk as scaled with its stan- 
dard deviation STk for GOE (a) and GUE (b). Events corresponding 
to Tk = were neither included in the distributions nor in the av- 
erage and the standard deviation used in the shifting and rescaling of 
Tk- The dashed line corresponds to a Gaussian with zero mean and 
unit variance. The solid line corresponds to Eq. for J/D = 0.24. 



Random Matrix Theory. For g ^> 1, the energy levels 
within an energy window smaller than E c obey the Wigner- 
Dyson distribution of RMT 1 13]. The eigenenergies s n and 
eigenfunction intensities x n fluctuate independently. The lat- 
ter are themselves uncorrelated and obey the Porter- Thomas 
distributions I3l 1 1 311 - Figure ^ shows the distribution of Tk 
obtained from the unfolded spectrum of random matrices of 
size N — 500 using 500 realizations for each Gaussian en- 
semble, namely orthogonal (GOE) and unitary (GUE). The 
Fermi energy is taken to be in the middle of the band. The de- 
pendence of the average (Tk) on J at (Tk) > A is found 
to coincide with the clean case result, which is given by 
T^ « 0.57 .D exp(— D/ J) for a spectrum of equally spaced 
levels. For small J, Tk approaches A; in a system without 
fluctuations of energy levels and wavefunctions, it would turn 
abruptly to zero at J* » D/[ln(2iV)-f 0.58] in a nonanalytical 
fashion, T$\j -► J*) = -A/In [{D/J* - D/J)/A] /2. 
However, fluctuations are important at (Tk) < A, leading 
to non-zero values for Tk below the clean-limit threshold 
J*. In Fig. \l\ the distributions of Kondo temperatures for 
the GOE and GUE are shown for different values of J. The 
Kondo temperatures are shifted and rescaled in order to facil- 
itate the comparison with a Gaussian distribution. The latter 
might be expected from a naive use of the central limit the- 
orem: Since wave function amplitudes corresponding to dis- 
tinct eigenstates are uncorrelated in RMT, Tk is determined 
by a sum of independent random variables when energy level 
fluctuations are neglected. However, even when (Tk) exceeds 
A there are deviations from the Gaussian behavior. This is 
in agreement with approximate analytical results obtained by 
performing the average over the wave function amplitudes ex- 
actly but assuming equally spaced energy levels and Tk ^> A, 
when we find in the limit of large number of levels N in a 



P(Tx)~-exp 

_ rp(0) 



[3k In 



feK \ 



(3n 



(2) 



where n — T K /A and (3 denotes the ensemble (j3 = 1 
for the GOE and 2 for the GUE). This distribution yields 
(T K ) = T ( °\ independent of (3. ¥oxT K « T { ° } , it is approxi- 



mately Gaussian with a width STk b = y A T^ / (3, in agree- 
ment with an earlier derivation 1 1 511 - but note that the distri- 
bution is not log-normal as was incorrectly concluded, there. 
The ratio between the variance in the orthogonal and unitary 
regime is thus \/2 for t£?' > A, in agreement with numerical 
results 1 14]. When time-reversal symmetry is present, the en- 
ergy level repulsion is weaker and the tendency to localization 
stronger, resulting in a wider distribution than in the unitary 
regime 1 16]. Within RMT, the width is found to vanish in the 
thermodynamic limit, when A — > 00. 
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FIG. 2: (Color on-line) The distribution of Tk with eigenstates of 
a 20 x 20 Anderson model without magnetic field (solid symbols) 
compared with the distribution for the GOE with N = 400 (dashed 
line). J/D = (a) 0.1, (b) 0.35, and (c) 0.6. For the Anderson model, 
4,000 realizations were used for each value of g. 

Anderson model. Going beyond RMT, stronger deviations 
from the Gaussian behavior can be expected due to correla- 
tions in the wave functions. In order to describe the conduc- 
tion electrons, we use the noninteracting tight-binding Ander- 
son model with nearest-neighbor hopping t and random site 
energies drawn independently from a box distribution of width 
W centered at zero. We assume each single-electron eigen- 
state to be spin degenerate and consider 20 x 20 square lattices 
with periodic boundary conditions. The Hamiltonian was di- 
agonalized and the eigenenergies and eigenvectors were used 
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to determine Tjc through Eq. Q without unfolding of energy 
levels. The disorder amplitude W and the parameter g are 
related in the Born approximation by 1/g = jgW 2 /(Ept), 
where the lattice constant is a = 1 (N = L 2 ) and A = 8t/N. 
This expression is valid for weak disorder g > 1 and for a 
fiat density of states. We chose Ep = 2t to avoid traces of 
the van Hove peak at the band middle. Only time-reversal 
symmetric cases where studied. The resulting distribution of 
Tk is shown in Fig. |2] where a non-Gaussian shape is evi- 
dent. The GOE and Anderson model distributions only agree 
at weak disorder and for J <C D, Fig. |2ja). For larger J 
[see Fig. |2jb)], the distribution remains non-Gaussian for all 
disorder strengths, while the GOE distribution corresponding 
to the same value of J/ D is almost Gaussian. The disordered 
Anderson model presents a bimodal structure in P(Tk) that is 
absent in the RMT case. That structure is strongly pronounced 
in the strong disorder regimes, with two distinct peaks clearly 
visible. As J/ D increases further, weight is transfered from 
the low-TV to the high-TV peak [see Fig. EJc)]. 

Analytical information on the fluctuations of Tjc can be ob- 
tained by relating it to the correlation function of the LDOS 
0. Defining by Eq. Q with average DOS v , we can 
rewrite Eq. Q in terms of the LDOS as [17] 
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Expanding TV to second order in 5p = p — v, we can re- 
late it to the correlation function of LDOS, i?2 DOS M = 
(p(r, E) p(r, E + uj)) /v 2 . While in RMT there are no cor- 
relations between wave functions at different energies, in the 
disordered Anderson model they are of order 1/g, 

L? (|^„(r)| 2 |^ m (r)| 2 ) ee f nm ( u ) = 1 + |Ren w , (4) 

independently of the states n,m (n ^ m). The correla- 
tion is stronger in the time-reversal symmetry case. The de- 
pendence on disorder is determined by a sum over diffuson 
modes, yielding, for L > I and in two-dimensions (2D), 
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In terms of the spectral correlation function, R 2 (uj) = 
(p(E)p(E + lo))/v 2 - 1, the correlation function of LDOS 
reads 



Rf os {io) = R 2 {lo) f nm (oj) + 6(u/A) f nn (cu). 



(6) 



For lo < E c , R2(u>) has an oscillatory correction of order 
1/g 2 to the leading RMT term |3j. For w > E c , the oscilla- 
tory part of the spectral correlation function decays exponen- 
tially, while there is a correction of order 1 / g 2 which decays 
like 1/lo without oscillations |3]. Expanding the exponent to 
second order in Sp, we find 
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where Sp arises from the spectral self-correlation term propor- 
tional to S(u)/A) and is of order A/Tk- To this order we also 
obtain (T K ) = T [ ° ] + {6T K ) 2 /2T i ° ) . The decaying parts of 
the spectral correlation function and the correlations of wave 
functions yield 
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In the thermodynamic limit E c — > 0, the leading term in Eq. 
(|8} arises due to the correlations between wave functions; to 
leading logarithmic order, we obtain 



C (t,T k ) 



1 



■In 



3 ( Ep 



6nPg \gT K 



(9) 



which is valid for D/J^S>\ng. Thus, the width of the distri- 
bution is finite in the thermodynamic limit and the deviation of 
the average Kondo temperature from the clean limit is stronger 
than expected from Anderson's theorem when applied to the 
Kondo problem, in contrast to an earlier prediction 1 18J . A 
finite STk will have implications to the thermodynamic prop- 
erties of disordered metals with magnetic impurities, as well 
as to heavy fermion materials. An earlier work 1 19] consid- 
ered the effect of fluctuations of LDOS as well, but obtained 
an even larger effect by disregarding that wave function cor- 
relations are only of order 1/g. 
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FIG. 3: (Color on-line) Dependence of average Kondo temperature 
(Tk) (a) and standard deviation (b) in the Anderson model on the di- 
mensionless conductance g for different values of J/D. Conditions 
are identical to Fig. |2| with averages taken over 500 realizations. 
The solid lines correspond to Eq. keeping E c and A finite. No- 
tice that since we assume E c t < 1, the lines are only presented for 
1< g < 20. 

The expressions for (Tk) and STk based on Eq. @ are 
compared with results from the numerical simulations in Fig. 
[5] They are found to agree qualitatively with the dependence 
on the disorder strength g seen in the simulations for different 
values of J. It is only at g > 1 and (T K ) > A that (T K ) 
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is insensitive to disorder. In general, for weak disorder, there 
is good agreement between GOE and Anderson model results 
for STx (not shown). At stronger disorder, fluctuations in the 
Kondo temperature increase in amplitude due to the increase 
of correlations in the LDOS at different energies over an in- 
terval of order of 1/r, Eq. (0}. 
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FIG. 4: (Color on-line) Probability of free moments as a function 
of J for the Anderson model at different disorder strengths g (same 
conditions of Fig. Also shown are the numerical result for the 
GOE with TV = 400. J* was estimated for a flat band. 

Free moments. The long tail toward small-T^ indicates 
that disorder enhances the probability of having free magnetic 
moments at low temperature. This is confirmed in Fig. 0] 
where we show the probability of finding no solution to Eq. 
([0 as a function of J/D. We note that for T K < A only 
the two closest levels to the Fermi energy are strongly cou- 
pled to the magnetic impurity. In the clean limit, the proba- 
bility of free moments is a sharp step function of J, namely, 
Pfroo(J < J*) = 1, and Pf ree (J > J*) = 0. When disorder 
is present, J* fluctuates due to the random wave function in- 
tensities x n and level spacings s n of the two closest levels to 
the Fermi energy 11411 . The probability that the level spacing 
at the Fermi energy is of the order of (J) is exponentially 
small for J > J*. However, in Fig. |4]we see that the decay 
with J is slower than that, indicating that the fluctuations of 
wave functions are crucial in the appreciable increase in the 
probability of free moments for J > J*. 

In the absence of an external magnetic field, magnetic im- 
purities can break time-reversal symmetry since the spin dy- 
namics of the magnetic impurities is slow compared to the 
time scale of the conduction electrons [20]. For a mesoscopic 
sample at T > E c , weak localization corrections are sup- 
pressed when A"W L = 1/E c t s exceeds 1 $2$. At finite tem- 
perature, the spin scattering rate l/r s is renormalized by the 
Kondo correlations: It is small at both T > T K and T <CT K , 
having a maximum at around T = Tk 12 ill , which is given by 
a fraction a of the unitary limit of the scattering cross sec- 
tion. In 2D, one finds l/r smax = ariMVF / '(2/cf), where 
riM is the concentration of magnetic impurities, with a as 0.2 
12 2l l24ll . Thus, the maximum value of the crossover param- 
eter is given by X s m = cxNm/9, where Nm = n^jL 2 is 
the number of magnetic impurities in the sample. When there 
are only a few magnetic impurities, Nm < g, X s < 1, the 



metallic grain is in the orthogonal regime. Increasing the con- 
centration of magnetic impurities increases the parameter X s , 
thus decreasing the width of the distribution of Kondo tem- 
peratures. The maximum in the temperature dependence of 
the spin scattering rate results in a plateau of the dephasing 
time lEHl . Mesoscopic Au wires with Fe impurities at con- 
centrations ns = 3—10 ppm have a Kondo temperature 
T K « 300 mK J2l|. While both E c and A are found to 
be of order 10~ 7 K, and thus completely negligible in such 
wires, we can estimate that J/D = 0.074 is much larger than 
J*/D = 0.035. Therefore, we obtain from Eq. (|9} that the 
width of the distribution of Kondo temperatures STk is of the 
same order as (T K ), with (T K ) = 300 mK. 

To conclude, the distribution of Kondo temperatures has 
a finite width in disordered metals due the correlations of 
eigenfunctions within an energy interval of width 1/r. Thus, 
the presence of free magnetic moments in disordered metals 
at temperatures below the average Kondo temperature is ex- 
pected to yield a finite contribution to the electron dephasing 
rate and to the thermodynamic properties of disordered metals 
in the presence of localized magnetic moments 1 23 ] . 
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Note: After the submission of this work we learned that a 
result similar to Eq. (|9jl was recently obtained independently 
by Micklitz and collaborators |25]. After the submission we 
learned of the work by Cornaglia and collaborators who also 
studied numerically the distribution of the Kondo temperature 
with the same approach, and found a peak at small Tk, as well 
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